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Abstract 

We generalize to three active flavours a previous two flavour model for the resonant 
spin flavour conversion of solar neutrinos to sterile ones, a mechanism which is added 
to the well known LMA one. The transition magnetic moments from the muon and 
tau neutrinos to the sterile play the dominant role in fixing the amount of active 
flavour suppression. We also show, through numerical integration of the evolution 
equations, that the data from all solar neutrino experiments except Borexino exhibit a 
clear preference for a sizable magnetic field either in the convection zone or in the core 
and radiation zone. This is possibly related to the fact that the data from the first 
set are average ones taken during a period of mostly intense solar activity, whereas 
in contrast Borexino data were taken during a period of quiet sun. We argue that 
the solar neutrino experiments are capable of tracing the possible modulation of the 
solar magnetic field. Those monitoring the high energy neutrinos, namely the fiux, 
appear to be sensitive to a field modulation either in the convection zone or in the core 
and radiation zone. Those monitoring the low energy fluxes will be sensitive to the 
second type of solar field profiles only. In this way Borexino alone may play an essential 
role, since it examines both energy sectors, although experimental redundancy from 
other experiments will be most important. 
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1 Introduction and Motivation 



Although the effort in solar neutrino investigation has decreased in recent years, several 
intriguing questions in this area remain open. Their clarification may lead to a better 
knowledge of the neutrino intrinsic properties, the structure of the inner solar magnetic 
field, or possibly both. In fact it is still unclear for example whether the active solar 
neutrino flux varies in time [1-4] or why the SuperKamiokande energy spectrum appears to 
be flat [5,6]. The generally acknowledged LMA solution [7] to the neutrino deficit observed 
by all solar neutrino experiments does not explain these facts, while it further predicts an 
event rate for the Chlorine experiment [8] which is 2a above the observed one [9]. These 
could be indications of physics beyond LMA. 

Based on this motivation and in line with the originally proposed resonant spin flavour 
precession of solar neutrinos [10,11], we were lead to develop a model whereby neutrinos are 
endowed with a transition magnetic moment converting active into sterile ones by virtue 
of its interaction with the magnetic field of the sun [12]. In particular we considered a 
scenario in which the apparently time varying event rate of the Gallium experiments [13, 14] 
was viewed in connection with the solar magnetic activity [15] . Owing to the uncertainties 
involved, such data are however also consistent with a constant Ga rate, which is the alter- 
native conventional view we will consider in the present paper. It is also possible to generate 
time variations of the active neutrino flux using the parametric resonance for matter den- 
sity perturbations in the presence of a radiation zone magnetic field without resorting to 
magnetic moments or sterile neutrinos [16]. 

The model expound in refs. [12, 15] considered two resonances, the LMA one between 
two oscillating active neutrinos and the spin flavour precession one determined by a tran- 
sition moment between one of the active flavours and the sterile one. The location of the 
active-sterile resonance, near the bottom of the solar convection zone, was fixed by the cor- 
responding active sterile mass squared difference which for this purpose was chosen to be 
O (10~^ — 10~^) eV^. The survival and transition probabihties were calculated using the 
Landau-Zener approximation. 

In this paper we extend our previous model with two active flavours and one sterile to 
the more reahstic case of three active flavours and a sterile. Whereas the Landau Zener 
approximation works well in the LMA resonance, this is not so for spin flavour precession, 
thus we will resort to the numerical integration of the evolution equations. However the best 
fit parameter values are obtained from the Landau Zener approximation, as they arc found 
to coincide in both methods. We take several values of ^^13 in the allowed range for both 
strong and weak solar fields. The model event rates for all solar neutrino experiments are 
evaluated and confronted with the data. Special emphasis is given to the SuperKamiokande 
energy spectrum [6] and the recent energy spectrum from the Borexino experiment [17]. 
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We consider two classes of solar field profiles, whose field strengths peak in the solar core 
(Wood-Saxon potential type) and in the lower convection zone. Together these span a large 
number of possibilities. 

The paper is organized as follows: in section 2 the derivation of the (4x4) Hamiltonian 
in the mass eigenstate basis is presented. In subsection 3.1 the already known fact that the 
survival probability is a decreasing function of ^13 is shown to follow from a simple argument. 
In subsection 3.2 the evaluation of the event rates and spectra is presented. The results are 
discussed in subsection 3.3. Finally in section 4 our conclusions are drawn. Our numerical 
calculations are based on the updated central values [18] for Am^^, ^12, ^23, Am|2 and we 
use a neutrino transition moment between flavour states not larger than /ii, — 1.4 x 10"^^//^. 
As for we chose to investigate three cases: ^13 = 0, 0.1 and the central value, 0.13. The 
fits to all data, including rates and spectra (except for Borexino) improve once the magnetic 
field is introduced. As regards Borexino, the fit worsens in this case. This may be connected 
to the fact that all former experiments report time averaged data taken during times of more 
or less intense solar magnetic activity, whereas Borexino data were taken during a period of 
minimum activity. In contrast, solar data alone show no clear preference for a vanishing or 
sizable 613. 



2 The Hamiltonian 



The (4x4) Hamiltonian involving one sterile neutrino and three active ones will be expressed 
in the mass eigenstate basis. This is related to the flavour basis by 
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where U^^^^{AxA) is the straightforward (4x4) extension of the usual leptonic mixing 
matrix [19]. As before [12, 15] no vacuum mixing between the sterile and any of the active 
flavours is assumed, so that the free propagating term of the Hamiltonian is in the mass 
basis 
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and the matter (interaction) term is in the flavour basis 
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Here ^[e,n,T)s are the transition magnetic moments between the active flavours and the sterile 
one and B is the magnetic flcld proflle. The quantities Vc-, Kt are the refraction indexes 
for charged and neutral currents, namely Vc = GpV^Ne, Ki = —{Gp/ V^)Nn with 
(electron density) and Nn (neutron density). Given eq.(l), the mass and flavour Hamiltonian 
representations are related by 
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so that the full Hamiltonian in the mass basis is 
f Eo \ / 
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Subtracting Ei and denoting by /ii,2,3 the transition magnetic moment between mass eigen- 
states and 1,2,3 respectively, we have 
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Straightforward matrix algebra leads now to 
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where Ug- denotes the first row entries of the (3x3) JJ^^^^ matrix. In the following we 
assume vanishing phases. Eq. (7) is the mass basis Hamiltonian that we use throughout. 
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3 Probability and Rates 



3.1 3 Flavour probability and ^13 

A simple argument shows that the electron neutrino survival probability for three active 
flavours [20,21] 

Psxsil^e Z^e) = COs' ^13i'2x2(z^e M^e ^ iVe COs' ^13)) + siu^ ^13 (8) 

where extra terms 0(10^^) or smaller have been neglected [21], is a decreasing function of 
^13. In eq.(8) the (2x2) probability with the replacement A^e cos^ ^13 is given by [21] 

^2x2 = I + ^ cos 2^,2 (9) 
^ ^ V^(%^sin2M^ + 02 

where x is the fractional solar radius and 

0=^Ar,cos2^i3-^cos2^i2. (10) 

Straightforward calculations show that the derivative of ^3x3(1^6 —* v^) with respect to 6'i3 
is negative for all solar neutrino energies if 

P2X2 - ^ COS 20,, cos^ ^,3 3 - tan^ ^13 > 0. (11) 

4V2 [(^sin2^i2)2 + 02]t 

In fig.l we plot the quantity on the left hand side of this inequality as a function of the 
neutrino production point denoted by its fractional radius x, for energies = 0.1, 8, 18.8 
MeV. Nf. is evaluated from the data on p and Xh given in [22] with 

K = ^i±^ (12) 

where p, rrip and Xh are the density, the proton mass and the hydrogen mass fraction. 
All neutrino parameters including ^13 were fixed at their best fit central values [18]. The 
minimum of the quantity (11) is clearly seen in fig.l for the three energy values considered 
and the central value ^13 = 0.13 claimed in ref. [18]. For x = t/Rq, this minimum reaches 
zero at a; ~ 0.2 as the energy is increased up to Emax = Emax{hep) = 18.8 MeV and 613 
up to 0.47 which is much above 3a from its central value 0. Thus we can conclude that for 
all experimentally allowed values of the physical quantities involved the condition (11) is 
satisfied and ^3x3(^^6 ^e) is a decreasing function of 6^13. 



-'^It has been recently pointed out [23,24] that the hint for a non-vanishing 6*13 is at most a la effect on 
the basis of present data. 
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3.2 Field profiles and rates 



We base our numerical calculations for the rates on the standard solar model with high 
metalicity, BPS08(GS) [25]. As far as the solar field is concerned, solar physics provides 
very limited knowledge on its magnitude and shape. For instance in ref. [26] upper bounds 
of 0.5-1.5 G and 30 G are quoted in the bottom of the convection zone and in the core 
respectively while in ref. [27] an argument is presented favouring an upper bound of 600 G 
in the radiation zone. On the other hand, the authors of ref. [28] estimate in the bottom 
of the convection zone an upper limit of 300 kG and in the mid-radiation zone and solar 
centre a magnetic field strength of 0.7 MG and 7 MG respectively [29]. 

Given the above mentioned uncertainties we consider the two following plausible profiles 
which are approximately complementary to each other (see also fig. 2) 

Profile 1 

B = — , ° < X < 0.71 (13) 

cosh[6(x - 0.71)] ^ ' 

B = ° nr^.M 0-71 < 2: < 1 (14) 

cosh[15(x - 0.71)] ^ ' 

Profile 2 

B=- r-^- -r 0<a;<l, (15) 

1 + exp[10(2x - 1)] ' ^ ^ 

Profile 1 has a peak Bq at the bottom of the convection zone, for fractional solar radius 
X ~ 0.71, its physical motivation being the large gradient of angular velocity over this 
range [30]. It should not exceed 300 kG at this depth and 20 kG at 4-5% depth, hence its 
fast decrease along the convection zone [28]. Profile 2 is of the Wood-Saxon type, being 
maximal at the solar centre. In this case the peak field Bq could be as large as a few MG [31]. 

The survival probability was obtained from the numerical integration of the Schrodinger 
like equation with Hamiltonian (7) using the Runge-Kutta method. All event rates (total 
and spectral) were evaluated as described in refs. [32,33]. The expression for the ^B spec- 
tral rate as applied to the SuperKamiokande and Borexino experiments with three active 
neutrinos is now 
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where CTg is the charged and neutral current cross section and o"^ is the neutral current one. 
The energy resolution functions are given in [34] and [35] and the threshold energies E^. — 5 
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MeV, Ef. = 2.8 MeV for SuperKamiokande and Borexino respectively. For the statistical 
analysis of all solar data (except Borexino) we used the standard definition [32, 33] 



where indices ji, j2 run over solar neutrino experiments and the error matrix includes the 
cross section, the astrophysical and the experimental uncertainties. 

The rates we obtained were confronted with the data from the CI experiment [8], the 
Ga ones [13], the SuperKamiokande spectrum [6], the SNO rates and spectra [36], the 
Borexino spectra for neutrinos [17] and for the remaining fluxes [37] Starting with 
profile 1 (fig.2) which peaks at the bottom of the convection zone, we considered the case 
of a relatively strong field and a vanishing one, and likewise for profile 2 (fig.2). Numerical 
results are insensitive to field values below 50 kG. In order to provide a feeling of the rates 
variation with 613, we also run this parameter from zero to 0.13. The values of the active 
sterile mass squared difference, which determines the location of the spin flavour resonance, 
were Amgi = ml — mf = 1.25 x 10"''' eV^ and Amg^ = 2.7 x 10~^ eV^ obtained by fltting 
with proflle 1 and proflle 2 respectively, using the Landau Zener approximation for the spin 
flavour precession resonance. For the flux normalization we used Jb = 0.95 [17] and the 
remaining neutrino parameters were [18] 



Am^i = 7.67 X 10-^eV^ Aml^ = 2.39 x 10" sin^iz = 0.559, sin^as = 0.683. (18) 



Since, as explained below, our predictions refer to the average magnetic activity in the 
solar cycle, they should be evaluated for an intermediate fleld strength. Therefore we take 
Bq = 140 kG for proflle 1, while for proflle 2 the peak value will be obtained from fltting. 
Regarding the mass basis magnetic moments /ii,2,3; we found that the numerical results 
are independent of which moment is chosen to be the largest. For the assumed peak fleld 
strength, the flts require this largest value to be 2 x To this end the two flavour 

(transition) moments IJI'{^,t)s must be of order 1.4 x 10"^^/is with equal or arbitrarily 
smaller. A value 1.7 x lO"^^//^ is also possible for the largest of the mass basis magnetic 
moment^, provided either /i^^ or /i^^ is as large as 1.4 x 10"^^/!^ with the other two no 
smaller than 1.0 x lO"^^/!^. In other words, /i(^,T)s are the dominant moments that flx 
the amount of the active neutrino flavour suppression. In the following we will consider 
fi(^^r)s — 1-4 X 10^^^/iB with fies arbitrarily equal or smaller, thus ensuring that the largest 
of the /i's is ^ 2 X lO^^^yU^- 

^We recall that, as mentioned in the introduction, we neglect any possible time variation in the Ga rate 
and perform fits to its average value, see ref. [13]. 
■^This would however require Bq ~ (160 — 170) kG. 




(17) 
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3.3 Discussion 



Starting with profile 1, the results for all solar neutrino experiments except Borexino, 
which is discussed below, are displayed in table I (total rates and fluxes) and flg.3 (Su- 
per Kamiokande spectral rate). They show a clear preference for a sizable magnetic field. 
We note that not only the fiatness of the spectrum is enhanced thus providing a better fit 
to the data (see fig. 3), but also the total rate predictions for the CI, SK and SNO experi- 
ments strongly improve (table I). A strictly constant spectrum could on the other hand be 
obtained by varying the solar parameters ^m\^ and sin ^12 within a 2a range. As for the 
Ga rate, vanishing and sizable fields are equivalent, as both classes of predictions lie within 
la of the central value [13]. Moreover as to the magnitude of ^^13, the predictions do not 
show any clear preference. Altogether the results are more sensitive to changes in the solar 
parameters than in the atmospheric ones. 
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64.3 
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1.70 


2.25 


0.95 


38.4 


44.1 


83.4 



Table I - Peak field values (profile 1), sinOi^, total rates (in SNU for Ga and CI experiments, 
in lO^cm^^s^^ for SK and SNO), and the corresponding 's- The total number of degrees 
of freedom is 82 = 84 experiments ( Ga + CI + 44 SK + 38 SNO data points) - 2 parameters, 
(see ref. [15]). It is seen that for a sizable field (Bq = 140 kG) all fits improve. 



The Borexino spectral rate for is shown in fig. 4. The top curve is the central value of 
the ratio between the best fit recoil spectrum due to oscillated neutrinos and the spectrum 
due to non-oscillated ones evaluated from fig. 3 of ref. [17]. The next group of three curves 
represents the predicted spectra for vanishing field with ^13 = 0, 0.1, 0.13 from top to bottom 

respectively, and the bottom three curves represent the same but for a field {Bq = 140 kG). 
The result of a analysis with 4 degrees of freedom for vanishing field and profile 1 with 
Bq = 140 kC is presented in the first three columns of table II. Owing to the magnitude of 
the experimental errors, it is not possible to conclude whether the data shows any preference 
for a vanishing or a finite 613, although a sizable magnetic field appears to be disfavoured. 
Prom fig.4 and table II our predictions look more sensitive to the magnetic field variation 
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than to the ^13 one within their respective allowed ranges. An improved significance can 
be obtained if Borexino are able to reduce their errors to 1/3 of the present ones (see the 
fourth column of table II) . 
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Table II - The result of our analysis for Borexino (profile 1): from the first three columns 
it is seen that no conclusion can he drawn as for the magnitude of sindi^ and that the 
significance is too low for the data to favour a vanishing field. The last column shows the 
variation if the experimental error were reduced to 1/3, so a vanishing field would clearly 
he favoured. The result would still in this case he inconclusive regarding the size of sinOi^ . 

Prom the data and the model predictions presented, it is unclear whether Borexino 
can favour either a neghgible or a sizable solar magnetic field. The available data from the 
remaining experiments appear to favour on the other hand a sizable field, possibly connected 
to a more intense solar activity in the convection zone. Hence it might be appropriate to 
examine the period in which the data were taken. In particular the SuperKamiokande 
spectrum refers to the period from May 31, 1996 to July 15, 2001 during which the average 
sunspot number was 65. On the other hand the Borexino spectrum refers to a data 
taking period from July 15, 2007 to June 21, 2008 with average sunspot number 4 [38]. In 
most of the former period the solar magnetic activity increased and reached an 11-year peak 
in the Summer of 2000, whereas in the latter the activity was continuously at its minimum. 
Therefore in the light of the present model, one expects the present Borexino spectrum for 

to coincide with the LMA prediction and the SuperKamiokande one to reflect an active 
sun. 

We have also tested the model for the remaining fluxes observed in the first Borexino 
phase [37], including in addition the pp and pep neutrinos. In this case £" < 1.7 MeV for all 
fluxes, so that with AmQ^^ = 1.25 x 10~^ eV^ all neutrino resonances lie below x = 0.5 where 
the magnetic field strength is -B < ^Bq (see fig. 2). Furthermore since the matter density 
is larger in this range, the variation of the field besides being smaller, becomes much less 
important. The maximum variation that is reflected in the event rates is no greater than 
1%, hence, given the order of magnitude of the experimental errors which is about 25% [37], 
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it cannot be expected to be seen in this case for any flux. We may therefore conclude that 
it is of prime importance that Borexino will continue monitoring both the low energy and 
the flux during the present increasing solar activity period. 

While it is generally accepted that the sunspot activity is interrelated with the possible 
modulation of the convection zone magnetic field, no connection appears to exist between 
such varying activity and the radiation and core magnetic field. There is however a recent 
claim in the literature [2] suggesting the existence of an inner tachocline separating the 
core from the radiation zone and an inner dynamo producing a strong magnetic field and a 
second solar cycle. Independently of the fact that our knowledge on this matter that can 
be obtained from solar physics is very limited at present, it will be shown in the following 
that the solar neutrino data are consistent with a varying field in these inner regions. 

Referring to profile 2 [eq.(15) and fig. 2] and in order that a possible time modulation 
may be detected, all resonances from active to sterile must be located deeper inside the 
sun with relation to profile 1, so that the shift from a weak to a strong field or vice versa 
is reflected in the intensity of the neutrino flux. In this way the best fit to the data was 
obtained for Am^Q = 2.7 x 10~^ eV^ and Bq = 0.75 MG □ with the remaining neutrino 
parameters as in eq.(18). The results for the total rates and fluxes are now shown in table 
III and those for the SuperKamiokande and ^B Borexino spectra in figs. 5 and 6. 
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Table III - 


Same as table I for 


profile 


2 where 


the vanishing field 


case is 


omitted. As for 



profile 1, with a sizable field (Bq = 0.75 MG) all fits improve with relation to the vanishing 
field (compare with table I). 

Again, as for profile 1, the SuperKamiokande data show a clear preference for a large 
field (fig. 5) and the quality of the fits is the same for both profiles. As can be seen from 
a comparison between figs. 3 and 5 or figs. 4 and 6, the spectra for a sizable peak field are 
much alike and it will be very hard to experimentally distinguish between them in this way. 
The actual difference can be explicitly seen in fig. 7 where the Borexino spectrum for both 
profiles is shown for 613 = with the remaining parameters as in tables I and III. Whereas 

^For instance with Am^Q = 1.25 x 10"^ eV^, the neutrinos with energy E = 5 MeV have their resonance 
at a; = 0.64 whereas for Am?n — 2.7 x 10^^ eV^ this resonance moves to x = 0.34. 
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for profile 1 the spectrum presents a shallow minimum around E — 8 MeV, it decreases 
monotonically in the case of profile 2. We have also seen that in the case of profile 2 the 

results are less stable, in the sense that small variations either in the neutrino parameters 
Am2i and 612 or Bo lead to larger variations in the rates and fits. Similarly to profile 1, the 
results for a analysis for profile 2 are shown in table IV and the case for a vanishing or 
a sizable field is again inconclusive. 
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Table IV - The same as table II for profile 2. The fourth column indicates the variations 
relative to its values for vanishing field given in the first three rows of table II. 

Finally for profile 2 we have calculated the rates corresponding to the remaining fiuxes 
which were observed in the first Borexino phase [37]. For all these neutrinos, with E < 1.7 
MeV, the corresponding resonances lie in the range x < 0.23 where the field is close to 
maximal. Thus the event rate modulation is expected to be much stronger than with 
profile 1. It is now approximately 9% which we believe to be within reach of the Borexino 
experiment in the future. The results are shown in fig. 8 where we chose to represent the 
modulation of the ^Be, ^^O and ^^N fiuxes. Owing to the magnitude of the errors involved 
(~ 25%), the of order 38 (51 d.o.f.), is extremely fiat for both profiles and it is hard to 
distinguish any preference of the data at all either for profile 1 or 2. So it is also of prime 
importance to keep Borexino monitoring the low energy neutrino fiuxes too, namely ^Be 
and CNO as in the first phase. Therefore solar neutrino experiments hold the potential to 
clearly trace a field modulation inside the sun and moreover possibly to distinguish whether 
this modulation occurs in the convection zone or deeper in the radiation zone and core. 



4 Summary and Conclusions 

We have developed a model with three active neutrino fiavours communicating to a sterile 
one in matter with magnetic field through transition magnetic moments. Its motivation is 
to provide better fits to the solar neutrino data than the LMA ones, in particular a flat 
SuperKamiokande spectrum and a better prediction for the Chlorine rate, while keeping 
accurate predictions for all other rates including the recent ^B Borexino spectrum. We in- 
vestigated two magnetic solar field profiles, one which peaks at the bottom of the convection 
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zone and another at the solar centre. These represent two classes of plausible possibihties 
which somehow complement each other. 

The starting point was the derivation of the (4x4) Hamiltonian followed by a simple 
and general argument showing that the survival probability is a decreasing function of the 
still unknown mixing angle ^13. This fact, which is reflected more strongly in the charged 
current data, leads to a parallel shift of the spectral event rates. This shift is however not 
enough to distinguish a clear preference of the data between a vanishing or a sizable ^13. 

We found that among the three transition moments, the ones connecting and 1/^ to 
the sterile are the dominant ones that fix the amount of active flavour suppression. They 
are both required to be of order 1.4 x IO^^^/ib, while z/gs may be equal or arbitrarily smaller. 
Alternatively either z/^ or fj^r separately may be of order 1.4 x 10~^^//b with the remaining 
two of order 1.0 x 10"^^//^ in which case a slightly stronger field is required. 

On the other hand it was found that all experimental data, with the exclusion of the 
Borexino ones, favour a relatively large magnetic field of either class. To this end it is 
important to realize that the former data are average ones and refer to extended periods. 
In particular the SuperKamiokande spectrum refers to a period when the average solar 
magnetic activity was relatively intense, and hence it is sensible to expect it to be flat in 
a way that it reflects a large field in accordance with the model predictions. As regards 
the Borexino spectrum in a similar energy range, it is not possible at present to conclude 
whether it favours the LMA spectrum or the LMA one with spin flavour precession, as the 
data errors are too large. 

Whereas the neutrino fluxes observed in the first Borexino phase were found to be 
insensitive to field modulations in the convection zone (profile 1), this is not so if the field is 
concentrated in the core and radiation zone (profile 2). The event rate modulation expected 
in this case is of the same magnitude as the one expected for the flux with any of 
the profiles. Hence we believe it extremely important to keep Borexino taking data for all 
neutrino fluxes during at least the first half of the present solar cycle expected to peak in 
2011 or 2012. 

Our claim is not that there is evidence of variability of the solar field proflle in the 
convection zone or cquivalcntly in the core and radiative zone, but rather that the neutrino 
data arc consistent with the possibility of either phenomenon. 

Our results concerning field profiles and data variability are qualitatively summarized in 
table V: a magnetic field concentrated around the bottom of the convection zone like profile 
1 can only show its modulation through an experiment monitoring the high energy fiux, 
whereas a field concentrated in the core and radiation zone like profile 2 can be detected by 
experiments monitoring either the high energy ^B or the low energy fiuxes. 

To conclude, solar neutrino experiments may hold a non-negligible potential to ascertain 
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whether there is a varying magnetic field inside the sun possibly connected to solar activity, 
a fact which otherwise may be very difficult to establish on the basis of solar physics alone. 
Moreover we have shown that it is also possible to trace whether this varying field is lying 
mostly at the bottom of the convection zone or deeper in the core and radiation zone. 



Varying field 


fiux 


Others 


Profile 1 (CZ) 


Yes 


No 


Profile 2 (WS) 


Yes 


Yes 



Table V - The possibility for detecting through solar neutrino experiments the magnetic 
fields concentrated either in the convection zone (profile 1) or in the core and radiation zone 
(profile 2). 
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Figure 1: The three lines represent expression (11) as a function of the neutrino production 
point for neutrino energies 0.1 MeV, 8 MeV and 18.8 MeV with sin6'i3 = 0.13. For all 
experimentally allowed values of 6*13 the quantity plotted in this graph is positive implying 
that the survival probability is always a decreasing function of 9i^. 
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Figure 2: The two solar field profiles given by eqs.(13), (I4) and (15) normalized to their 
peak field values and expressed as a function of the fractional solar radius. 
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Figure 3: The SuperKamiokande spectrum: theoretical predictions and data points [6] nor- 
malized to BPS08(GS) [25]. The top three curves refer to sin ^13 = 0,0.1,0.13 from top 
to bottom in the case of zero magnetic field, and the lower three curves refer to the same 
values of sinOis for a sizable field (profile 1), with B = 140 kG at the peak. There is a clear 
preference for a sizable field possibly related to solar activity, in comparison to a vanishing 
one. Units for the observed electron energy E^. are in eV . 
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Figure 4: The Borexino spectrum normalized to BPS08(GS) [25]. The top curve is 
extracted from fig. 3 of ref. [17] as explained in the main text. The lower two groups are from 
top to bottom the model predictions withsmOi^ = 0,0.1,0.13, B = and sin ^13 = 0,0.1,0.13, 
B = 140 kG at the peak (profile 1). There is a preference in this case for a vanishing field 
possibly related to a quiet sun. Notice that the theoretical curves coincide for > 5 MeV 
with the corresponding ones in the previous figure apart from a minor difference due to the 
energy resolution functions. 
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Figure 5: Same as fig. 3 with the three bottom curves referring to profile 2 with B = 0.75 
MG at the peak. 



20 



1.2 
1.0 
0.8 
0.6 



0.4 - 



0.2 - 



0.0 



-| 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 

0.47 I I I I I I I I I I 



0.37 







,9.6x10^ I.OxlO^H 



1 



J I I I I I I I I I I I I I I I I I L 



3.0x10^ 5.0x10^ 7.0x10® 9.0x10® 1.1x10'' 

E, (eV) 



Figure 6: Same as fig.4 (^B Borexino spectra) for profile 2. 
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Figure 7: Borexino spectra for neutrinos evaluated for profiles 1 and 2 at the best fit 
with = (parameter values as in the main text). The spectrum for profile 1 exhibits a 
shallow minimum while for profile 2 it is monotonically and smoothly decreasing with the 
energy. 
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Figure 8: Borexino spectra for'^ Be neutrinos (full lines), ^^O (dashed) and ^^N (dot-dashed) 
evaluated for vanishing field and profile 2 at the best fit with = (parameter values as 
in the main text). 
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